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Abstract 

The total cross sections for the elastic electroproduction of £ and J/^/j mesons for 
Q 2 > 8 GeV 2 and (W) ~ 90 GeV/c 2 are measured at HERA with the HI detector. 
The measurements are for an integrated electron— proton luminosity of ~ 3 pb _1 . 
The dependences of the total virtual photon— proton (j*p) cross sections on Q 2 , W 
and the momentum transfer squared to the proton (t), and, for the p, the dependence 
on the polar decay angle (cos#*), are presented. The J ftp : g cross section ratio is 
determined. The results are discussed in the light of theoretical models and of the 
interplay of hard and soft physics processes. 
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1 Introduction 



The study of elastic production of vector mesons in photo- and leptoproduction (Fig. |l|a) 
in fixed target experiments has provided information on the hadronic component of the 
photon and on the nature of diffraction. 

With the advent of the electron-proton collider HERA, there is renewed interest in 
vector meson production, in particular at large Q 2 (Q 2 is minus the square of the ex- 
changed photon four- momentum). HERA experiments have observed in deep-inelastic 
scattering that the proton structure function F 2 increases rapidly |l] with increasing W, 
the 7*p invariant mass, which is in striking contrast with the slow rise of the total 7p 
cross section at Q 2 ~ [§[. In addition, the rise of F 2 already at relatively low Q 2 values 
(Q 2 $ 2 GeV 2 ) indicates that the transition between these two behaviours is rapid ||. The 
study of the elastic production of vector mesons is expected to provide useful information 
about these different regimes, and in particular about the transition between them. 

e(k') e(lc') 




Figure 1: Diffractive vector meson production: a) elastic; b) proton dissociation. 

The subject of this paper |4| is the study of elastic g and J/ip meson electroproduction 
at large Q 2 (Q 2 > 8 GeV 2 ) and high W ((W) ~ 90 GeV/c 2 ), in the reactions 

e p — > e g p; g — > n + tt~ , (1) 

e p — > e J/ip p; J/ip —* Z + / ~; I = e or /i. (2) 

After the presentation of the event selection and of the mass distributions, the total 
cross sections and the differential distributions for reactions (|]) and (0) are studied. 
Features relevant for the study of the transition regime are then discussed, in particular 
the energy dependence of the cross section, the t distribution (t is the square of the four- 
momentum transferred from the photon to the target), the vector meson polarisation, 
and the evolution of the J/ip '■ Q cross section ratio. The results are discussed in the light 
of the interplay of hard and soft physics processes. Results on g production in a similar 
kinematic range have been presented by the ZEUS Collaboration 0. 

2 Models and phenomenology 

The elastic production of light vector mesons, in particular of g mesons, by real or quasi- 
real photons (Q 2 ~ 0), a process called hereafter photoproduction, exhibits numerous 
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features typical of soft, hadron-like interactions. These are a small angle peak in the 
distribution of the vector meson scattering angle with respect to the incident photon 
beam direction ("elastic peak"), a steepening of this distribution with increasing energy 
("shrinkage"), for W £ 10 GeV/c 2 a slow increase with energy of the cross section, and 
s-channel helicity conservation (SCHC). 

These observations support the vector meson dominance model (VDM) [/J, accord- 
ing to which a photon with energy greater than a few GeV behaves predominantly as 
the superposition of the lightest (q, w, 0) J PC = I— mesons. In this framework, the 
cross section for elastic vector meson production is related to the total meson— proton 
cross section through the optical theorem, and the energy dependence is related to the 
"universal" energy dependence of the total hadron— proton cross section. Neglecting a 
contribution ("reggeon exchange") which decreases with energy approximately as s~ ' 5 , 
the latter is parameterised at high energy as a tot oc s s with 5 ~ 0.0808, y/s being the 
hadron— proton centre of mass energy ||. This is expressed in Regge theory as due to 
soft pomeron exchange. The elastic production of g mesons has been studied extensively 
by fixed target experiments in photoproduction 0, [9| , for intermediate Q 2 |1(J and for Q 2 
£ 6 GeV 2 [O. Also at HERA, the photoproduction of g mesons exhibits the characteristic 
features of soft interactions [0 . 

In contrast, the VDM approach does not give a satisfactory description of J /if) pho- 
toproduction data |13|, [L4|, pL5|| : the cross section is smaller than the VDM prediction (see 
IH) and it increases significantly faster with energy than expected from soft pomeron 
exchange. This is also observed by the HERA experiments [ lTfl . 

Using the hard scale provided by the mass of the charm quark, the fast increase of the 
J /if) photoproduction cross section was predicted in the framework of QCD by Ryskin 
||18|| . In his model, the interaction between the proton and the cc pair in the photon 
is mediated by a gluon ladder, and non-perturbative effects are included in the nucleon 
parton distribution. 

At high energy and high Q 2 , g meson production is modelled using two different 
approaches based on QCD, which both refer to the pomeron as basically a two gluon 



system (see [|j|]). The g meson production is thus related to the gluon distribution in the 



proton, but the two approaches emphasize respectively soft or hard behaviour. 



In the soft model initially proposed by Donnachie and Landshoff fl20| , pi] , 16], the gluons 
are described non-perturbatively and the elastic cross section is expected to increase slowly 
with energy: da/dt (t = 0) oc W AS , with 5 ~ 0.0808. The vector mesons are predicted to 
be mostly longitudinally polarised, the cross section to fall as Q~ 6 and the J /if) cross 
section to be comparable to that of the g. 

In the second approach, perturbative QCD calculations similar to the work of Ryskin 
have been performed by several groups [22], ^3], |24fl , a hard scale being provided by the 
photon virtuality. A major prediction of this model is a rapid increase of the cross sec- 
tion with W, as a consequence of the rise of the gluon distribution in the proton. It is 
stressed, however, that non-perturbative processes are also expected to be present in an 
intermediate energy region |2^, ^] . The scattering amplitude is obtained in these calcula- 
tions from the convolution of the hadronic wave function of the photon, of the scattering 
amplitude of this hadronic component, and of the final state vector meson wave function. 
This is because the photon is viewed at sufficient energy as the coherent superposition of 
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hadronic states formed well before the target (essentially qq pairs), whereas the final state 
meson is formed beyond it, the corresponding time scales being much longer than the 
interaction time. As a consequence, the spatial dimensions of the hadronic wave function 
are an essential parameter in the interplay of perturbative and non-perturbative effects. 
It is predicted that the hard effects should show up earlier for small size objects than 
for large ones, and for longitudinally than for transversely polarised photons. The Q 2 
dependence of the j*p cross section is predicted, as in the non-perturbative approach, to 
fall as Q~ 6 but when the evolution of the parton distribution and quark Fermi motion are 
taken into account, it was pointed out that the Q 2 distribution is expected to be harder 
2|. In view of the more compact J/i[) wave function, the J/ip : g cross section ratio has 



been predicted |25| to exceed, at very high energy, the value 8 : 9 obtained from SU(4) 



and the quark counting rule |2C 



In addition to these two "microscopic" approaches inspired by QCD, calculations are 
also performed for elastic vector meson production on the basis of multiple pomeron 
exchange, with the effective pomeron intercept depending on the photon virtuality J27]] . A 
stronger increase of the cross section with energy is again predicted than for soft pomeron 
exchange with 5 ~ 0.0808. 



3 Detector and event selection 

The data presented here correspond to an integrated luminosity of 2.8 pb _1 for the g and 
3.1 pb _1 for the J/ip mesons. They were collected in 1994 using the HI detector. HERA 
was operated with 27.5 GeV positrons and 820 GeV protonsQ. The detector is described 



in detail in ref. 28 



The event final state corresponding to reactions (|IJ) and (Q) consists of the scattered 
positron and two particles of opposite charges, originating from a vertex situated in the 
nominal e + p interaction region. In most cases, the scattered proton remains inside the 
beam pipe because of the small momentum transfer to the target in elastic interactions. 

In the Q 2 range studied here, the positron is identified as an electromagnetic clus- 
ter with an energy larger than 12 GeV, reconstructed in the backward electromagnetic 
calorimeter (BEMC) and associated with a hit in the proportional chamber (BPC) placed 
in front of it at 141 cm from the nominal interaction vertexF]. The BEMC covers the 
polar angles 151° < 9 < 176°. Its electromagnetic energy resolution is cte/E ~ 6 — 7% 
in the energy range of the positrons selected for the present studies. The BPC angular 
acceptance is 155.5° < 6 < 174.5°. The scattered positron polar angle 9 e is determined 
from the positions of the BPC hit and of the interaction vertex. The trigger used for 
the present analyses requires the presence of a total energy larger than 10 GeV deposited 
in the BEMC, outside a square of 32 x 32 cm 2 around the beam pipe. Additional cuts, 
similar to those used for the structure function analysis [la], are applied to the hit and 
cluster position and shape in order to provide high trigger efficiency and good quality 

1 For the g studies only positron data are used, whereas the small amount of data taken with e~p 
scattering is included in the J/ip sample; in this report, positrons refers both to positrons and electrons. 

2 The forward (+z) direction, with respect to which polar angles are measured, is defined as that of 
the incident proton beam, the backward direction is that of the positron beam. 
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positron measurement. These cuts are complemented by the selection of events with Q 2 
> 8 GeV 2 . 

The decay pions (g events) or leptons (J /if) events) are detected in the central tracking 
detector, consisting mainly of two coaxial cylindrical drift chambers, 2.2 m long and 
respectively of 0.5 and 1 m outer radius. The charged particle momentum component 
transverse to the beam direction is measured in these chambers by the track curvature 
in the 1.15 T magnetic field generated by the superconducting solenoid which surrounds 
the inner detector, with the field lines directed along the beam axis. Two polygonal 
drift chambers with wires perpendicular to the beam direction, located respectively at 
the inner radius of the two coaxial chambers, are used for a precise measurement of 
the particle polar angle. For the present analyses, two tracks with transverse momenta pt 
larger than 0.1 GeV/c are required to be reconstructed in the central region of the tracker. 
For q production, the polar angles must lie in the range 25° < 9 < 155°, corresponding 
to particles completely crossing the inner cylindrical drift chamber for interactions at 
the nominal vertex position. For J/ip production, the accepted range is extended to 
20° < 9 < 160° in order to increase the statistics as much as possible while keeping good 
detection efficiency. The vertex position is reconstructed using these tracks. No other 
track linked to the interaction vertex is allowed in the tracking detector, except possibly 
the positron track. To suppress beam-gas interactions, the vertex must be reconstructed 
within 30 cm of the nominal interaction point in z, which corresponds to 3 times the 
width of the vertex distribution. In addition, the accepted events are restricted to the 
range 40 < W < 140 GeV/c 2 for g mesons and 30 < W < 150 GeV/c 2 for J/ip mesons. 

The tracking detector is surrounded by a liquid argon calorimeter situated inside the 
solenoid and covering the polar angular range 4° < 9 < 153° with full azimuthal coverage. 
In the case of elastic interactions, the calorimeters should register only activity associated 
with the decay particles or the positron. However, due to noise in the calorimeters and to 
the small pile-up from different events, elastic g and J / ip production can be accompanied 
by the presence of additional energy clusters. These are considered in terms of the variable 
Emaxi defined as the energy of the most energetic cluster which is not associated with a 
track. The E max distribution shows a peak at small values, attributed mostly to elastic 
and diffractive interactions, and a broad maximum at higher energies. The cut E max < 1 
GeV is applied to enhance exclusive g production. This is discussed in section || together 
with the effect of the cut \t\< 0.5 GeV 2 . 

The sample of events containing a positron and a vector meson candidate includes 
two main contributions: elastic production (see Fig. [TJa), defined by reactions (JIJ) and 
(0), and events where the proton is diffractively excited into a system X of mass Mx, 
which subsequently dissociates (Fig. |l|b). Non-resonant background is also present. It 
is possible to identify most of the "proton dissociation" events with the components of 
the HI detector placed in the forward region p9| , namely the forward part of the liquid 
argon calorimeter 4° < 9 < 10°, the forward muon detectors (arrays of muon chambers 
placed around the beam pipe in the proton direction, 3° < 9 < 17°) and the proton 
remnant tagger (an array of scintillators placed 24 m downstream of the interaction point, 
0.06° < 9 < 0.17°). When particles from the diffractively excited system interact in the 
beam pipe and the collimators, the interaction products can be detected in these forward 
detectors. The events are tagged as due to proton dissociation by the presence of a 
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cluster with energy E^ r larger than 1 GeV (0.75 GeV for the J /if) candidates) at an 
angle 9j^ r < 10° in the liquid argon calorimeter, or by at least 2 pairs of hits in the 
forward muon detectors (one hit pair is compatible with noise), or by at least one hit in 
the proton remnant tagger. This last criterion is not used for the J /if) events, since they 
have a flatter t distribution than the g events and would be partially vetoed by the proton 
tagger. 

In order to minimise the effects of QED radiation in the initial state, the difference 
between the total energy and the total longitudinal momentum E — p z of the positron 
and the two particles emitted in the central part of the detector is required to be larger 
than 45 GeV. If no particle, in particular a radiated photon, has escaped detection in the 
backward direction, E — p z should be twice the incident positron energy, i.e. 55 GeV. 

The selection criteria for the two samples, supplemented by the mass selections dis- 
cussed in section |5], are summarised in Table |1|. 





Q 


J/iP 


positron selection 


em. cluster > 10 GeV in 


BEMC outside 32 x 32 cm 2 




reconstructed positron energy > 12 GeV 




associated BPC hit 




Q 2 > 


8 GeV 2 


VM reconstruction 


2 tracks fitting to vertex (+ possibly e + ) 
p t > 0.1 GeV/c 




25° < 9 track < 155° 


20° < 9 track < 160° 




40 < W < 140 GeV/c 2 


30 < W < 150 GeV/c 2 


background suppression 


l^vertex Z 


nom. | *^ 30 Cm 




E max < 1 GeV 






\t\ < 0.5 GeV 2 




p. dissoc. ev. tagging 


E L A r < 1 GeV for 9 L A r < 10° 


E L A r < 0.75 GeV for 9 L A r < 10° 




< 1 hit pair in forward muon detectors 




no hit in proton tagger 




radiative corrections 


E- Pz 


> 45 GeV 


mass selection 


0.6 < m^ +n - < 1.0 GeV/c 2 


2.8 < m l+l - < 3.4 GeV/c 2 


selected sample 


180 events 


31 events 



Table 1: Selection criteria for g and J/ if) events. 



4 Kinematics and cross section definitions 

The kinematics of reactions (|1|) and (fj) are described with the variables commonly used 
for deep-inelastic interactions. In addition to s (the square of the e + p centre of mass 
energy), Q 2 and W, it is useful to define the two Bjorken variables y = p ■ q/p ■ k (in the 
proton rest frame, the energy fraction transferred from the positron to the hadrons) and 
x = Q 2 /2p ■ q, where k, p, q are, respectively, the four-momenta of the incident positron, 
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of the incident proton and of the virtual photon. These variables^ obey the relations 
Q 2 = xys and W 2 = Q 2 (± - 1). 

The kinematical variables can be reconstructed from four measured quantities: the 
energies and the polar angles of the scattered positron and of the vector meson. With 
the "double angle" method |30| used for the present analyses, Q 2 and y are computed 
using the polar angles 9 and 7 of the positron and of the vector meson, which are well 
measured: 

Q gmTjl + cosg) 

sm 7 + sin 6 — sin(7 + u) 
sin 6 (1 — cos 7) 
sin 7 + sin 9 — sin(7 + 0) ' 

where E is the energy of the incident positron. 

The meson momentum components are obtained from the measured decay products. 
The momentum of the scattered positron is computed from Q 2 and y, which provides 
better precision than the direct measurement. The energy transfer to the proton being 
negligible, the absolute value of t is given by: 

1*1 - (Ptp) 2 = ifte +Ptv) 2 , (5) 

where p tp , Pt e and p tv are, respectively, the momentum components transverse to the beam 
direction of the final state proton, positron and vector meson[]. 

The fourth quantity which is directly measured, the positron energy, is used to compute 
the variable E — p z : 

E~p z = (E e + E v ) - (p ze + p zv ), (6) 

E e and E v being the energies of the scattered positron and of the vector meson, and p ze 
and p zv their momentum components parallel to the beam direction. 

The cross section for elastic electroproduction of a vector meson V can be converted 
into a 7*p cross section using the relation 

d 2 (T t ot{ep -» eVp) , I/W1 , D v 

— ^ = T <Ttat{l P^Vp) = T a T (7 P -> Vp) (1 + e R), (7) 

where a to t, &T and ai are the total, transverse and longitudinal 7*p cross sections, 

R =cr L /a T , (8) 

and T is the flux of transverse virtual photons given by 

a em ( l-y + y 2 /2) 
Q 

e is the polarisation parameter 



Y _ "em ^ - y -r y 1 ^} _ ^ 



y + y 2 /2' 



(10) 



3 In this paper, the positron and proton masses are neglected. 

4 The lowest \t\ value kinematically allowed, t m i n ~ (Q 2 + my) 2 m 2 /y 2 s 2 , is negligible in this 
experiment. 
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Information on the vector meson production process can be obtained from the angular 
distributions of the decay particles. In particular, the probability for the g meson to 
be longitudinally polarised can be determined from the distribution of cos#*, where 6* is 
the angle, in the g rest frame, between the direction of the positively charged decay pion 
and the g direction in the 7*p centre of mass system (helicity frame) [[/], 31]: 



<1A ocl-r{ft + (3r$-l) cos 2 r. (11) 



dcos 9* 

With the assumption of s-channel helicity conservation (SCHC), is related to R: 



1 r 04 
1 '00 



E 1 — Tn 



oo 



5 Mass distributions and final samples 

Fig. shows for the selected events (Table [l|) the distribution of the invariant mass 
m^+Tr- in the range m n + n - < 2 GeV/c 2 , obtained by assigning the pion mass to the 
particles detected in the central tracker. The mass distribution without the E max and t 
cuts (insert in Fig. ga) is seen to peak at small m^^- values. The E max cut strongly 
reduces the background of events containing neutral particles, and enhances the g peak. 
The t cut is very effective in rejecting non-resonant events containing, in addition to the 
g candidate and the positron, particle(s) with a significant transverse momentum, which 
is not used to compute p tp in eq. (|J). This cut also enhances the elastic production 
signal compared to the background of proton dissociation events, which are known to 
have a flatter t distribution. In total, 180 events are found in the g peak region with 
0.6 < m w+w - < 1.0 GeV/c 2 . 

In Fig. ||a, events compatible with the (f> mass, when the charged particles detected in 
the central tracker are considered as kaons, have been removed {ttik+k- < 1-04 GeV/c 2 ). 
Assuming that vector mesons are produced according to the quark counting rule with 
the SU(3) ratios g : u : = 9: 1 : 2, the Monte Carlo simulation described below 
indicates that the remaining (p and u reflections contribute, in the m n + n - range (0.4 — 
0.6) GeV/c 2 , 2.0% of the g signal in the peak, and 0.7% in the range (0.6 - 1.0) GeV/c 2 . 
These contributions were subtracted statistically. 

The m 7T + 7T - distribution is described by a relativistic Breit-Wigner function, over a 
non-resonant background attributed to incompletely reconstructed diffractive photon dis- 
sociation. The Breit-Wigner function has the form |H| 

(13) 



dN(m 7T7T ) m nn m p r(m„ 



dm n7T (m 2 - m 2 J 2 + m 2 r 2 (m 7r7r ) ' 
with r(m 7r7r ) the mass dependent width 

r(m^) = T p A 3 2 - - . (14) 

Here m p is the g resonance mass and V p the width; q* is the pion momentum in the 
(7r + 7r~) rest frame, and q$ this momentum when m n + n - = m p . 
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Figure 2: a) 



mass distribution for the events with Q 2 > 8 GeV 2 , 40 < W < 140 



GeV/c , |i| < 0.5 GeV , E max < 1 GeV and no proton dissociation signal in the forward 
detectors; the superimposed curve is the result of a fit to a relativistic Breit-Wigner 
distribution over the background (15), which is described by the dashed curve. In the 
insert: same distribution without the E max and t cuts; b) mi+i- mass distribution for the 
events with Q 2 > 8 GeV 2 , 30 < W < 150 GeV/c 2 and no proton dissociation signal in the 
forward detectors; the superimposed curve is the result of a fit to a Gaussian distribution 
over an exponential background in the range 2.4 < mi+i- < 4.5 GeV/c 2 . 



The background has been parameterised using the distribution 

= «! (m„-2m T ) a2 e - Q3m ™, (15) 



djTTl/'jY'jY 

where is the pion mass and a±, «2 and are free parameters. This form, which 
includes a two pion threshold and an exponential fall off, is in qualitative agreement with 
the background shape in the insert of Fig. 0a. 

With these parameterisations, the resonance mass is 763 ± 10 MeV/c 2 and the width 
is 176 ± 23 MeV/c 2 , in agreement with the Particle Data Group (PDG) values of 770 
and 151 MeV/c 2 p3| . No skewing to low values of m^^- is needed to describe the g 
shape (for the Ross-Stodolsky parameterisation [34], the skewing exponent is found to be 
n = 0.3 ±0.5). 

Two alternative forms have been used for the resonance width: 

IW) = r p (^) 3 ^ (16) 



and 



r(m_) = r p (^) 3 . (17) 
% 



n 



They also give a g meson mass and width compatible with the PDG data. 

The non-resonant background under the g peak is estimated to be 11 ± 6%. The 
error includes the uncertainty on the resonance parameterisation and on the background 
shape, estimated by using for the latter an alternative linearly decreasing form. 

The distribution of the invariant mass rrii+i- for the selected events in the J / ip region 
is presented in Fig. ^|b. The J /ip mass is 3.13 ± 0.03 GeV/c 2 . The peak width is slightly 
larger than, but compatible with, the expectation obtained from the detector simulation. 
No t selection is applied since the t distribution is significantly flatter than for the g 
mesons (see Fig. §). No E max cut is required in view of the small background in this high 
mass region (compare Fig. |2|b and insert in Fig. |2|a): all J/ip candidate events have E max 

< 1.35 GeV, of which 4 have E max larger than 1 GeV. 

A sample of 31 J/ip candidate events is thus selected with the cut |m ; +/- — m^\ 

< 300 MeV/c 2 , where is the J/ip meson mass. The non-resonant background is 
estimated by fitting the sidebands using an exponential distribution and amounts to 
roughly 20% (6.8 events). No lepton identification is required, but 10 of the J /ip candidate 
events contain two identified electrons and 7 contain two identified muons. 

One event with E max > 1 GeV is a ip' — > J/ipir n candidate, with one identified muon 
and neutral clusters detected in the electromagnetic part of the liquid argon calorimeter, 
attributed to the interaction of photons from ir° meson decay. The measured J/ip mass 
is 3.16 ± 0.04 GeV/c 2 . The invariant mass computed using the two charged tracks and 
the neutral clusters is 3.67 ± 0.09 GeV/c 2 , in excellent agreement with the Particle Data 
Group value (3.69 GeV/c 2 ) |3|. 

Kinematical characteristics of the selected events are summarised in Table |2|. 





Q 




<Q 2 > [GeV 2 ] 


13.4 ±0.4 


17.7 ±1.5 


(W) [GeV/c 2 ] 


81 ±2 


92 ±6 


(scattered e + energy) [GeV] 


25.5 ±0.1 


24.9 ±0.3 


(scattered e + p t ) [GeV/c] 


3.4 ±0.1 


3.9 ±0.2 


(meson energy) [GeV] 


4.2 ±0.1 


6.0 ±0.3 


(tracks) [GeV/c] 


1.7 ±0.1 


2.3 ±0.2 



Table 2: Averages of kinematical variables characterising the selected events. 



6 Corrections and simulations 

Table ^| summarises the correction factors applied to the selected samples to take account 
of detector acceptance and efficiencies, smearing effects, losses due to the selection criteria 
and remaining backgrounds. 

Most corrections are estimated using a Monte Carlo simulation based on the vector 
meson dominance model, which permits variation of the Q 2 , W and t dependences, as 
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well as of the value of R P5| . The HI detector response is simulated in detail, and the 



events are subjected to the same reconstruction and analysis chain as the data. 

The accuracies of the g {J/ip) variable measurements are for W, 3.3 (4.3) GeV/c 2 , for 
Q 2 , 0.4 (0.4) GeV 2 , for t, 0.06 (0.10) GeV 2 . The values of the t slopes are little affected 
by the detector resolution. 

The scattered positron selection criteria induce Q 2 dependent losses for Q 2 < 12 GeV 2 . 
The error quoted in Table |3] corresponds to a systematic uncertainty on the positron 
direction of 2 mrad. The charged track selection criteria induce W dependent losses 
for small and high W, depending on the accepted W range in the two selections. A 
Q 2 — W correlation of the losses is observed, and taken into account in the corrections. 
The correction for the t cut in the g sample is computed using the measured t slope (see 
section 17.21) . 





Q 




trigger 

positron acceptance (Q 2 dep.) 

BPC hit — cluster link 

tracker acceptance (W dep.) 

track reconstr. (per track) 

track p tmin (per track) 

\t\ cut 

E — p z cut 

Efnax Cut 

forward det. cuts (\t\ dep.) 
mass selection 


1.01 = 
1.16 ±0.03 

1.03 = 
1.07 ±0.01 
1.05 ±0.03 

1.02 ±0.01 

1.03 ±0.02 

1.02 ±0.01 

1.03 ±0.03 

1.04 ±0.02 
1.22 ±0.01 


= 0.02 
1.15 ±0.03 

= 0.02 
1.29 ±0.03 
1.03 ±0.03 

1.00 ±0.01 

1.01 ±0.01 

1.03 ±0.03 
1.00 ±0.02 


non-resonant background 
proton dissoc. background 
cf) and uj background 


0.89 ±0.06 
0.91 ±0.08 
0.99 ±0.01 


0.78 ±0.14 
0.75 ±0.11 


photon flux / bin integration 
radiative corrections 
luminosity 


1.00 ±0.04 
0.96 ±0.03 
1.00 = 


1.00 ±0.07 
1.00 ±0.04 
= 0.02 



Table 3: Correction factors and systematic errors, averaged over the data samples. 

The choice of the E max cut for the g sample is a compromise between the loss of 
elastic events to which a cluster is accidentally associated in the calorimeter, and the 
presence of non-resonant background in the final sample. The loss is estimated using a 
Monte Carlo simulation which includes random noise in the calorimeters superimposed 
on elastic events. 

A simulation indicates that 2% of the elastic g events with |£] < 0.5 GeV 2 are lost 
because the proton has acquired sufficient p t to interact in the beam pipe walls, giving 
interaction products which are registered in the proton tagger; this loss is thus t dependent. 
For the J/ip sample, 1.5% of the events are lost because of interaction products registered 
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in the forward muon detectors. Another 2% loss for both samples is due to spurious hits 
in the latter. 

In view of the uncertainty on the high mass shape of the resonance, the g cross section 
is quoted in this paper for m n + w - < 1.5 GeV/c 2 , i.e. ~ m p + 5 Y p . The uncertainty is 
larger in the present case than for low energy data, for which a natural cut off is imposed 
by the limited available energy. With this definition of g resonance production line shape, 
the correction for the mass selection 0.6 < m n + n - < 1.0 GeV/c 2 is respectively 21% and 
23% for parameterisations (|T^) and (|TjJ). 

A simulation was performed in order to estimate the contribution to the final samples of 
proton dissociation events which are not tagged by the forward detectors. The distribution 
of the target mass Mx is parameterised as 1/M X . High mass states are assumed to decay 
according to the Lund string model |36| or, alternatively, to a final state with particle 
multiplicity following the KNO scaling law and isotropic phase space distribution. In 
the resonance domain, the mass distribution follows measurements from p dissociation on 
deuterium ||37| and resonance decays are described according to their known branching 
ratios. The decay particles are followed through the beam pipe walls and the forward 
detectors. 

For the g sample, the correction factor for the contamination of undetected proton 
dissociation events in the selected g sample is 0.91 ± 0.08. This number is obtained from 
the number of measured events tagged and not tagged by the forward detectors, and from 
the detection probabilities provided by the Monte Carlo simulation. No assumption needs 
to be made for the ratio of proton dissociation to elastic events. The error is a conservative 
estimate taking into account the uncertainties on the efficiencies of the forward detectors 
for tagging proton dissociation events and on the dissociation model. 

The correction factor for unobserved proton dissociation background in the selected 
J/ip sample, for which the proton tagger is not used, is 0.75 ± 0.11. 

The cross section measurements are given in the QED Born approximation for electron 
interactions. The effects of higher order processes are estimated using the HERACLES 



4.4 generator 38 



Radiative corrections for the g sample are of the order of 4% after the cut E — p z > 
45 GeV, and are weakly dependent on Q 2 and W. A systematic error of 3% is obtained 
by varying the effective Q 2 dependence of the 7*p cross section from Q~ A to Q~ 6 and by 
modifying the W dependence from a constant to a linearly increasing formQ. The small 
value of the correction is due to the high E — p z cut resulting from the good BEMC 
resolution; for the chosen value of the cut, small smearing effects are observed. 

For the J/ip sample, the radiative corrections determined using the measured Q 2 and 
W dependences of the cross section vary from +2% to —2%. 



7 Results 



5 In practice, the input to the program is an effective U F% structure function" parameterisation, with 
the chosen Q 2 and W dependences of the ep cross section for vector meson production. 
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7.1 Electroproduction cross sections 

The g and J /if) data are grouped in several {Q 2 , W) bins. Table [| gives, for each bin, 
the number of events, the integrated ep cross section and the 7*p cross section obtained 
using relation (171) for a given (Qq, VFo) value, taking into account the observed dependence 
across the bin. All known smearing, acceptance and background effects are corrected for. 
For the g sample, each event is weighted using the differential flux factor given by eq. (§). 
A 4% systematic error accounts for the uncertainty in the Q 2 and W dependences of the 
cross section used for the bin size integration and the bin centre correction. For the J/ip 
sample, in view of the small statistics, the photon flux is integrated over each {Q 2 ,W) 
bin. Since the data span a large range in Q 2 and W, this leads to a systematic error on 
the cross section of the order of 7%. 

The integrated cross section for g meson electroproduction with m^+^- < 1.5 GeV/c 2 

is 

a{e p -»• e g p) = 96 ± 7 (stat.) ± 13 {syst.) pb, (18) 

for Q 2 > 8 GeV 2 and 40 < W < 140 GeV/c 2 . 

The cross section for J/ip meson electroproduction, taking into account the J/ip — > 2 
leptons branching fraction 0.12 ||33|| , is 

o{ep-*e J/ip p) = 100 ± 20 {stat.) ± 20 {syst.) pb, (19) 

for Q 2 > 8 GeV 2 and 30 < W < 150 GeV/c 2 . 

The ratio J/ip : g is 0.64 ± 0.13 for Q 2 = 10 GeV 2 and 1.3 ± 0.5 for Q 2 = 20 GeV 2 . 

7.2 Momentum transfer distributions 

The t distributions of the selected g and J / ip events (Fig. |]) show the forward exponential 
peaking oc e bt characteristic of elastic interactions. 

For the g sample, the slope b of the t distribution is computed taking into account 
the contributions of the non-resonant and proton dissociation backgrounds estimated in 
sections |^ and Their exponential slopes were taken to be respectively 0.15 ± 0.10 
and 2.5 ± 1.0 GeV -2 , which is consistent with the t dependence of event samples which 
approximate these contributions. The results quoted below are rather insensitive to the 
choice of these slopes. The t distribution is corrected for detector effects, including the 
loss of elastic events tagged by the forward detectors. 

A fit for \t\ < 0.6 GeV 2 gives for the elastic slope the value b = 7.0 ± 0.8 ± 0.4 ± 0.5 
GeV" 2 (x 2 = 9.6 /10 d.o.f.) for Q 2 > 8 GeV 2 and 40 < W < 140 GeV/c 2 . The first error 
corresponds to the statistical precision of the fit. The second describes the spread of the 
fits according to the choice of the \t\ range (0.4 to 0.6 GeV 2 ) and of the E max cut value. 
The third comes from the uncertainties in the sizes and shapes of the backgrounds; it is 
dominated by the error on the total contribution of the non-resonant background. The t 
slope measured by the ZEUS Collaboration is b = 5.1 iJJ ±1.0 GeV~ 2 ||. The values of 
the slopes for two Q 2 and two W domains are given in Table [5]. 

For the J/ip sample, the slope value for \t\ < 1.0 GeV 2 is b = 3.8 ± 1.2 tfs GeV" 2 , 
after subtraction of the proton dissociation and non-resonant backgrounds with slopes 
b = 2 GeV" 2 . The systematic error is estimated by varying the background contributions 
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e p ^ e g p 




8 < Q 2 < 
40 < W < 80 GeV/c 2 


; 12 GeV" 
80 < W < 140 GeV/c 2 


number of events 
total correction factor 
integrated ep cross section [pb] 
Ql [GeV 2 ],^ [GeV/c 2 ] 
•y*p cross section [nb] 


57 

1.53 ±0.22 
31.5 ±4.2 ±4.5 
10, 65 

25.8 ±3.4 ±3.7 


47 

1.76 ±0.25 
29.8 ±4.3 ±4.2 
10, 115 

29.4 ±4.3 ±4.2 




12 < Q 2 ■ 
40 < W < 80 GeV/c 2 


< 50 GeV" 
80 < W < 140 GeV/c 2 


number of events 
total correction factor 
integrated ep cross section [pb] 
Q 2 [GeV 2 ], W [GeV/c 2 ] 
7*p cross section [nb] 


39 

1.32 ±0.19 
18.5 ±3.0 ±2.7 
20, 65 

5.0 ±0.8 ±0.7 


37 

1.23 ±0.18 
16.3 ±2.7 ±2.4 
20, 115 

5.0 ±0.8 ±0.7 


e p e J/ip p 




8 < Q 2 < 
30 < W < 90 GeV/c 2 


; 40 GeV 2 
90 < W < 150 GeV/c 2 


number of events 
total correction factor 
integrated ep cross section [pb] 
Ql [GeV 2 ], W [GeV/c 2 ] 
7*p cross section [nb] 


15 

1.50 ±0.27 
61 ± 18 ± 12 
16, 65 

7.8 ± 2.2 ± 1.6 


16 

0.92 ±0.17 
40 ± 12 ± 8 
16, 115 

12.2 ±3.4 ±2.5 




30 - lb < 
8 < Q 2 < 12 GeV 2 


150 GeV/c 2 
12 < Q 2 < 40 GeV 2 


number of events 
total correction factor 
integrated ep cross section [pb] 
Ql [GeV 2 ], Wo [GeV/c 2 ] 
7*p cross section [nb] 


10 

1.82 ±0.33 
49 ± 18 ± 10 
10, 88 

17.6 ±6.3 ±3.7 


21 

0.90 ±0.16 
51 ± 13 ± 10 

20, 88 

6.6 ±1.6 ±1.4 



Table 4: Numbers of events and cross sections for different Q 2 and W ranges, for g and 
for J/t/j elastic production; the 7*p cross sections are given for Q 2 = Ql, W = W . 
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Figure 3: \t\ distributions for a) the g sample; the data points are not corrected for 
background; the dashed line is the result of a fit taking the background into account, as 
described in the text; b) the J/ip sample, corrected for the presence of background. Both 
distributions are corrected for acceptance, losses and smearing effects. The solid lines 
correspond to the elastic exponential slopes. 



HI < 11' < 1 

8 < Q 2 < 12 GeV 2 


40 GeV/c 2 

12 < Q 2 < 50 GeV 2 


b = 7.8 ± 1.0 ±0.7 GeV~ 2 


b = 5.7 ± 1.3 ±0.7 GeV-' z 


8 < Q 2 < 
40 < W < 80 GeV/c 2 


50 GeV 2 

80 < W < 140 GeV/c 2 


b = 6.2 ± 1.0 ±0.7 GeV -2 


b = 8.0 ± 1.3 ±0.7 GeV-' z 



Table 5: Slopes of the g meson t distributions for different Q 2 and W domains. 
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by one standard deviation and their slopes between and 3 GeV -2 . The combined value 
of the HERA experiments [17a-b] for the slope in photoproduction is b = 4.0 ±1.0 GeV -2 . 
Three events have \t\ > 1.1 GeV 2 (see Fig. |3|b), of which one has E max > 1 GeV. These 
3 events contribute 15% to the cross section quoted in this paper. The ip' candidate 
event (see section |5|) has \t\ = 0.16 GeV 2 , t being computed including the neutral clusters 
attributed to tt° mesons (Q 2 = 26.3 GeV 2 , W = 72.6 GeV/c 2 ). 
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Figure 4: t slope for elastic g production in fixed target and HERA experiments a) as a 
function of Q 2 ; b) as a function of W for photoproduction and for Q 2 ~ 10 GeV 2 . 



Fig. |4]a shows that for the W domain of the HERA experiments the decrease with 
rising Q 2 of the t slope for g elastic production is similar to that observed at lower W. 

The comparison of the NMC and HI results for Q 2 ~ 10 GeV 2 (Fig. ^b) shows an 
increase of the t slope with energy. This shrinkage of the elastic peak with W (or y/s) is 



observed in diffractive hadron interactions [S7J and in photoproduction (see the compari- 



son of fixed target and HERA results in Fig. |]b). In the framework of Regge theory, for 
pomeron exchange and in terms of the exponential parameterisation, the shrinkage of the 
elastic peak can be written 

b(W 2 ) = b(W 2 = W/" 2 ) +2 a \n(W 2 /Wo), (20) 

where a' is the slope of the effective pomeron Regge trajectory: 

a P (t) = a P (0) + a' t. (21) 

Applying relation (EDI) to the Q 2 = 10 GeV 2 results (with statistical and systematic 
errors combined quadratically) gives for a' the value 0.41 ± 0.18 GeV~ 2 , in agreement 
with a value of 0.25 GeV -2 deduced from hadronic interactions For the HI data 

alone, there is also an indication for an increase of the slope with W (see Table |5|). 
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7.3 Q 2 dependence of the cross sections 

The Q 2 dependence of the total 7*p cross section (a tot = <r T + e a L ) for the elastic g meson 
production by virtual photons (Fig. |5]a) can be described by Q~ 2n with n = 2.5 ±0.5 ±0.2. 
In extracting the dependence on Q 2 of the cross section, correction has been made for 
the presence of non-resonant background, for which n = 1.5 ± 0.2 as obtained from 
the events with \t\ > 0.5 GeV 2 or E max > 1 GeV. The second error on n reflects the 
uncertainty on the background size and shape and the spread of the results according to 
the details of the fitting procedure. The cross section dependence for the present data 
is close to that obtained by NMC (n = 2.02 ± 0.07) and by ZEUS (n = 2.1 ± 0.4 t° i 
for 0.0014 < x < 0.004). It should be noted, however, that the NMC data span a large 
range in the polarisation parameter e from e = 0.50 at Q 2 = 2.5 GeV 2 to e — 0.80 for 
Q 2 > 10 GeV 2 , whereas the HERA data are for e = 0.99. Although the Q 2 dependence 
is probably sensitive to this kinematical effect, it was not taken into account because the 
evolution of R with Q 2 in the NMC data (see eq. |7]) is not published. The differences in 
the absolute normalisations are discussed in section [7.4j . 




10 20 30 



Q 2 [GeV 2 ] Q 2 [GeV 2 ] 

Figure 5: Q 2 dependence of the j^p — > Vp cross section a) for g production, the ZEUS 
points being given for the restricted x range 0.0014 < x < 0.004; b) for J/ip production 6 , 
the curves being the result of the fits described in the text. 

The Q 2 dependence of the J/ip production cross section at HERA is shown in Fig. [5|b. 
The errors on the high Q 2 data points include the uncertainty in the Q 2 dependence of 
the background. The evolution from photoproduction to high Q 2 is well described by 
1 / (Q 2 + m 2 b ) n with n = 1.9 ±0.3 (stat.). This is similar to the Q 2 dependence of the low 
energy EMC results [0, for which a fit of the data shown in Fig. [5]b gives n = 1.7 ± 0.1. 

7.4 W dependence of the cross sections 

The W dependence of the g (for m n + w - < 1.5 GeV/c 2 ) and J/ip production cross sections 
is shown on Fig. |6] for Q 2 = 10 and 20 GeV 2 . These values are chosen in order to minimise 
the bin centre corrections for the g analysis. 
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Figure 6: W dependence of the 7*-*- 1 p — > cross section in fixed target and HERA 

experiments a) for p production (computed for m n + n - < 1.5 GeV/c 2 ); b) for J/^ produc- 
tion. For the g data, an overall normalisation uncertainty of 31% for ZEUS and of 20% 
for NMC is not included in the plot. 
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The ZEUS results in Fig. have been scaled to the Q 2 values of the HI measurements 
using the dependence of the latter (see section |OD . These results can be directly compared 
to the HI results although they include no explicit cut-off on the m ir + ir - mass. Indeed, the 
cross sections quoted in || are determined assuming a non-relativistic Breit-Wigner mass 
distribution with mass independent width, integrated over the full kinematical range. It 
turns out that this procedure leads to a cross section closely similar to that obtained using 
the relativistic form of eq. (0-14) for m^+ v - < 1.5 GeV/c 2 , as is done in the present 



experiment. The ZEUS results are higher than those of HI, but the discrepancy is not 
very significant when the overall 31% systematic error on the ZEUS results ||, which is 
not included in the plot, is taken into account. Other differences between the results of 
the two experiments are observed: both the ZEUS t (section |7. 2|) and cos#* distributions 
(section |7.5|) are flatter than for HI. It is worth emphasizing in this context that the event 
by event selection in the HI analysis, which uses the forward detectors, provides a very 
clean elastic g sample. 

The NMC results shown in Fig. |6]a are obtained from the published Q 2 dependence of 
the cross section for interactions on deuterium (Fig. 4 of [lid]), which has small nuclear 
corrections POL O] . They are corrected for the different values of the polarisation param- 



eter e in the two experiments (the NMC measurement of R is used). The published NMC 
cross section was computed using for the g resonance the Breit-Wigner parameterisation 
given by eq. ( |13"D and (|16|), for m^+ n - < 1.5 GeV/c 2 (40j. The results of the two experi- 



ments can thus be directly compared. An overall 20% systematic uncertainty [lid] is not 
included in the plot. 

A significant increase with energy of the elastic 7*p cross section is observed from the 
NMC to the HERA domains. Following section 0, it is parameterised as da/dt (t = 0) oc 
W AS . The fitted values of 8 are for g elastic production: 

Q 2 = 10GeV 2 : 5 = 0.14 ±0.05, (22) 

Q 2 = 20GeV 2 : 5 = 0.10 ±0.06. (23) 

The errors result from the combination of statistical and systematic errors of both exper- 
iments, including the 20% normalisation uncertainty for NMC. 

The measurements ( p2|) and (^) take into account the following effects: 

— the da/dt (t = 0) cross sections are obtained by multiplying the total cross sections 
by the corresponding b slopes. The HI slopes given in Table |5| were used and the cor- 
responding NMC slopes were computed^ according to the shrinkage description given by 
eq. with a' = 0.25 GeV~ 2 . 

— the cross section definitionf] contains a kinematical factor due to phase space in- 
tegration, involving the centre of mass energy W, the mass squared of the particles and 
Q 2 . This factor is not part of the study of the interaction dynamics contained in the W 
evolution of the matrix element. As the Q 2 values considered here are rather large com- 
pared to W 2 for the NMC experiment and small for HI, there is a rising contribution to 
the W dependence amounting to 12% (26%) between NMC and HI energies for Q 2 = 10 

e For Q 2 = 10 GeV 2 , the use of the measured NMC slope instead would lead to an additional increase 
of S by 0.03. No measurement of the NMC slope is published for Q 2 = 20 GeV 2 . 
7 See eq. (23.32) and (23.36), p. 1292 of p§. 
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(20) GeV 2 . This corresponds to a decrease of S by 0.02 (0.03), which is included in the 
measurements ([22] ) and (|2~3|). 

Additional effects may have to be taken into account: 

— model predictions are often computed for the longitudinal cross section and not 
for the total cross section <jt + s a^. Taking into account the difference in the R values 
for the two experiments (see section |7.5| ), the 5 values for o L alone would be increased 
by 0.02 with respect to the values ( ^2|) and (^) (a possible Q 2 dependence of R is not 
considered) . 

— at the NMC energies, reggeon exchange could contribute significantly to elastic g 
production. Following the parameterisation obtained by Donnachie and Landshoff for 
the forward amplitude (eq. (9) of |16[) and assuming it holds for high Q 2 , the contri- 



butions to the da/dt (t = 0) cross sections of the purely reggeon exchange and of the 
reggeon— pomeron interference term are, respectively, 4% and 28% of the pomeron ex- 
change contribution (0% and 4% at HERA). To extract the forward differential cross 
sections from the measured total cross sections, assumptions have to be made concerning 
the relevant b slopes. For the purely reggeon exchange term, the value 0.83 is used for 
the slope of the Regge trajectory |fE? |, with the same parameter 6(W 2 ) = 2.5 GeV~ 2 as 



for the pomeron (Wo = 1 GeV/c 2 ). The slope for the interference term is chosen as the 
average of the pomeron and the reggeon slopes. When these contributions are subtracted, 
the value of S is increased by 0.02. 

Not including an error for theoretical uncertainties, the values of S for ol and pomeron 
exchange only are 0.18 ± 0.05 for Q 2 = 10 GeV 2 and 0.14 ± 0.06 for Q 2 = 20 GeV 2 . 

The W dependence of J /if) production is presented in Fig. [|bf] for Q 2 ~ 0, 10 and 20 
GeV 2 . A steep increase of the photoproduction cross section is observed from low energy 
to the HERA experiments. For higher Q 2 values, a similar increase is observed between 
the EMC and the HI measurements. However, quantitative comparisons should be taken 
with caution in view of normalisation uncertainties and the possible presence of inelastic 
background in the fixed target data. 

7.5 q decay angular distribution 

The acceptance corrected cos 9* distribution for the selected g sample is shown in Fig. f?]a. 
After subtraction of the non-resonant background, which is consistent with being flat in 



cos#*, and correction for detector effects, the fit of eq. (PD to this distribution gives 



,.,oi 

00 



0.73 ± 0.05 ± 0.02. The first error is statistical, the second reflects the uncertainty on the 
background subtraction. Assuming SCHC, relation ( |l"2"|) gives R = o~l /ot = 2.7 l^s -al> 
with (e) = 0.99. 

The value of R is shown in Fig. 0b together with fixed target measurements^ and 
given in Table |6| for two values of Q 2 and of W . Compared with results at low Q 2 , a clear 
increase of R with Q 2 is observed. 



An attempt was made to test the hypothesis [24| that the g meson should be completely 



longitudinally polarised for \t\ ^> Aq CD , by dividing the data with \t\ < 0.5 GeV into 

8 All results presented in Fig. ||b have been rescaled to take into account the latest measurements of 
the J/tp branching fractions: B(J/tp -> e+e") = 5.99 ± 0.25%, B(J/tp -> m + aO = 5.97 ± 0.25% @. 

9 EMC measurements [11a] with largest errors have been omitted. 
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Table 6: Values of R — &l /<jt for g elastic production in different Q 2 and W domains. 
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two samples with \t\ respectively smaller and larger than 0.15 GeV 2 . The predicted effect 
was not observed, but the imposed cut is probably too low to provide a sensitive test of 
the prediction. 



8 Discussion and conclusions 

The production of elastic g and J/i/j mesons by virtual photons has been measured at 
HERA with the HI detector. Samples of 180 and 31 events, respectively, have been 
collected with Q 2 > 8 GeV 2 and 40 < W < 140 GeV/c 2 (30 - 150 GeV/c 2 for the 
J/ip), for an integrated luminosity of 2.8 (3.1) pb -1 . Most of the proton dissociation 
background is removed using the forward components of the HI detector and the small 
residual backgrounds are corrected for. 

A major interest of the study of elastic vector meson production is that predictions 
have been proposed for the differential cross sections both in the framework of a soft, non- 
perturbative approach, and based on perturbative QCD calculations for hard processes. 

The main difference between the predictions of the soft and of the hard approaches 
concerns the rise of the cross section with energy, expected respectively to be slow or 
fast. For g production at high Q 2 , the W dependence of the cross section attributed 
to pomeron exchange is parameterised as da/dt (t = 0) oc W 45 . Using the NMC and 
the present HI results, the measured values of 8 for the total cross section or + £ &l 
are 8 = 0.14 ± 0.05 for Q 2 = 10 GeV 2 and 8 = 0.10 ± 0.06 for Q 2 = 20 GeV 2 . For the 
longitudinal cross section alone, and taking into account possible reggeon exchange, 
these values would be 8 = 0.18 ± 0.05 and 8 = 0.14 ± 0.06, respectively. 

For the soft pomeron model of Donnachie-Landshoff, the expected value is 0.08. For 
the hard approach, it is presumably in the range 0.20 — 0.25 [3b]. The present measure- 
ments thus lie between the values expected for these two types of models. 

The suggestion P31 that the g cross section measurement may provide information on 



the gluon distribution in the proton is applicable only when the hard regime is reached. As 
this condition does not seem to be fulfilled for the present W and Q 2 ranges, an attempt 
to extract the gluon distribution from these data seems premature. 

In contrast, the J/ip production cross section for Q 2 > 8 GeV 2 increases strongly 
from the fixed target to the HERA region. This increase is of the same order as in 
photoproduction. This indicates that a hard regime is reached for J/ijj production already 
at low Q 2 , which could be related to the smaller spatial extent of the wave function and 
the large scale provided by the charm quark mass. 

A major result of the present measurement is the similarity of the cross sections for 
g and J /if) elastic production. Whereas J /if) photoproduction, which is suppressed by 
factors of 100 to 1000 with respect to g, is not well described by the "quark counting 
rule", quark flavour symmetry appears to be approximately restored for Q 2 of 10 to 20 
GeV 2 . Such a behaviour is expected both in the soft and the hard models. However, 
this evolution is observed to be faster than for some hard models (a ratio 1/2 has been 
proposed for Q 2 ~ 100 GeV 2 |5|). 

The t dependence of the production differential cross sections for g and J/ip mesons 
are found to be well described at low t values by exponential dependences e bt . For the 
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g sample, b is 7.0 ± 0.8 (stat.) ± 0.4 (syst.) ± 0.5 (bg.) GeV -2 . This value is smaller 
than for photoproduction at HERA, showing that the decrease of the slope with rising Q 2 
observed in fixed target experiments extends to the HERA regime. This behaviour can 
be attributed to the decrease of the qq transverse separation in the photon with rising Q 2 . 

The evolution of the slope with W is sensitive to the interplay of soft and hard effects 
in q production. There is an indication that the shrinkage of the elastic peak observed in 
hadron interactions and g photoproduction also occurs in the present electroproduction 
data, at large Q 2 . This is as expected from Regge predictions based on soft pomeron 
exchange, in contrast with the little shrinkage predicted in perturbative calculations for 
hard processes |25| . 

For J /if) production, the t distribution is well described with b = 3.8 ± 1.2 j^jj GeV -2 , 
which is smaller than for the g. This difference can be qualitatively explained by the fact 
that the J /if} wave function is more compact than the g wave function. 

The Q 2 dependence of the g total cross section can be described by a power law Q~ 2n , 
with n = 2.5 ± 0.5 ± 0.2. This distribution is slightly harder than initially expected both 
for non-perturbative two gluon exchange and for hard QCD calculations (oc Q~ 6 )- It 
is compatible with predictions taking into account the Q 2 evolution of parton densities 
and the transverse motion of the quarks in the photon |25j. It should be noted that 



the calculations are performed for ol whereas the present measurement is of otou which 
includes a transverse contribution at the level of 20—25%, with presumably a steeper Q 2 
dependence p2|, p3| . 

The Q 2 dependence of the J/if> production cross section is parameterised as l/(Q 2 + 
rn 2 b ) n 1 with n = 1.9 ± 0.3 (stat.), which is similar to fixed target results. 

The polar angle distribution of the decay pions indicates that g mesons are mostly lon- 
gitudinally polarised: the spin density matrix element Tqq is 0.73 ± 0.05 ± 0.02. Assuming 
SCHC, R = a L /a T is 2.7 

-as -a!- The increase of R with Q 2 , which is predicted both 
by the non-perturbative model and the QCD calculations, is observed. The indication in 
the present data of an increase of R with W suggests, in the framework of hard physics, 
that perturbative features at high W and moderate Q 2 are indeed more important for 
longitudinal than for transverse photons. 

In conclusion, the study of g and J /if) meson production offers a contrasting picture. 

For J/if> mesons, hard physics effects are probably at work already for very small 
Q 2 . The present high Q 2 data support this interpretation, albeit with limited statistical 
precision. 

For g mesons, the W, t, Q 2 and cos 9* dependences of the cross section allow a more 
detailed study, from which a mixed picture emerges. It can be speculated that the present 
data correspond to a transition regime, with interplay of hard processes, amenable to per- 
turbative description, and soft processes, requiring a non-perturbative approach. The W 
dependence of the cross section does not provide conclusive evidence in favour of a purely 
soft or a purely hard model. The indication of shrinkage of the elastic peak with W is 
important because it shows a continuation with increasing Q 2 of the photoproduction 
behaviour and is at variance with expectations for a purely hard behaviour. The obser- 
vation that SU(4) flavour symmetry is restored in the J /if) : g cross section ratio, which 
is a striking feature of the present measurements, is expected in both models. The study 
of the Q 2 and cos#* distributions also does not discriminate between the soft and hard 
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models, since their predictions are similar. 
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